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In several classes of particle physics models - ranging from the classical Majoron models, to 
the more recent scenarios of late neutrino masses or Mass- Varying Neutrinos - one or more of the 
neutrinos are postulated to couple to a new light scalar field. As a result of this coupling, neutrinos 
in the early universe instead of streaming freely could form a self-coupled fluid, with potentially 
observable signatures in the Cosmic Microwave Background and the large scale structure of the 
universe. We re-examine the constraints on this scenario from the presently available cosmological 
data and investigate the sensitivity expected from the Planck satellite. In the first case, we find 
that the sensitivity strongly depends on which piece of data is used. The SDSS Main sample data, 
combined with WMAP and other data, disfavors the scenario of three coupled neutrinos at about the 
3.5a confidence level, but also favors a high number of freely streaming neutrinos, with the best fit 
at 5.2. If the matter power spectrum is instead taken from the SDSS Large Red Galaxy sample, best 
fit point has 2.5 freely streaming neutrinos, but the scenario with three coupled neutrinos becomes 
allowed at < 2a. In contrast, Planck alone will exclude even a single self-coupled neutrino at the 
4.2a confidence level, and will determine the total radiation at CMB epoch to AN^-^^ =io 3 i^'^ 
errors) . We investigate the robustness of this result with respect to the details of Planck's detector. 
This sensitivity to neutrino free-streaming implies that Planck will be capable of probing a large 
region of the Mass- Varying Neutrino parameter space. Planck may also be sensitive to a scale of 
neutrino mass generation as high as 1 TeV. 



I. INTRODUCTION 

Cosmology and particle physics are becoming increas- 
ingly linked, as is particularly evident in the dark sector, 
comprising dark matter, dark energy and dark radiation 
(neutrinos or other thermalized relativistic particles). 
Considerable work has been done to discover the nature 
of the new particle physics underlying the cosmological 
dark sector. While much remains unknown about it, the 
current abundance of data, combined with more precise 
measurements in the future, promise to give a clearer 
window into the dark. In this paper we consider the 
neutrino fraction of the dark sector. Cosmology already 
tells us a number of interesting facts about the neutrino 
sector when the universe was in its infancy. Big Bang 
Nucleosynthesis (BBN), for example, places significant 
constraint on the amount of radiation at T 1 MeV, 
translating to an effective constraint on the total num- 
ber of thermalized neutrino species, N;;ff < 4.5 The 
strongest upper bound on neutrino mass is also currently 
derived from the matter power spectrum on small scales, 
and the current limit is J2 < ^■^'^ at 95% C.L., as 
given in 0. 

From a particle physics perspective, there are 
many well-motivated models which generate beyond-the- 
standard- model neutrino physics, potentially giving rise 
to signals in the early universe. Most of these proposals 
could be classified into one of the following four broad 



categories: suggestions that neutrinos have (i) small but 
nonzero masses [1, 13] ; (ii) interactions with electromag- 
netic fields (via magnetic/transition moments) [5, _6]; (iii) 
interactions with new heavy states Q; (iv) interactions 
with new light states As is well known, the first of 

these ideas has actually proven true, giving us the first 
direct discovery of physics beyond the standard model. 
While it is certainly possible that there are no other ob- 
servable effects of new physics, it is intriguing to think 
that the neutrino sector could hold more surprises. The 
pragmatic approach is to set aside any theoretical prej- 
udice (which has proven rather unhelpful in the past in 
predicting neutrino properties) and settle the issue ex- 
perimentally. Should new neutrino interactions be dis- 
covered, the theoretical implications would be profound. 

Testing neutrino properties in terrestrial experiments 
has been very challenging. It took three decades from the 
first indication of the "solar neutrino problem" [l^] to 
conclusively demonstrate that neutrinos are massive and 
oscillate. As for the other proposals on the list, progress 
has been even slower. In the case of non-standard neu- 
trino interactions mediated by novel heavy particles, 
bounds from accelerator-based experiments solar 
[T3I [isll. a tmospheric [TgI . [17| and beam neutrino ex- 
periments [iSl . [l9| remain quite weak to this day. Bounds 
on the electromagnetic moments changed little over the 
last 15-20 years, with the best constraints coming from 
cooling of astrophysical systems such as red giants be- 
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fore helium flash [20|. Finally, while certain models of 
neutrino coupling to light scalar fields (those involving 
weak doublet [2111211 or triplet Q fields) have been ruled 
out by the LEP data on the Z boson width, a large class 
of models (invol ving coupling to singlet fields) remains 
viable (see, e.^., [23l|). 

In this paper, we show that the experimental sensi- 
tivity to nonstandard interactions mediated by a new 
light particle is about to be significantly improved. In 
this case the experimental tools are cosmological observa- 
tions rather than neutrino oscillation experiments. The 
new interaction can manifest itself through large scat- 
tering rates of the neutrinos with themselves (and the 
new light fields), or through extra radiation in the early 
universe (through thermally populating the new light 
states). As we will show, while the present cosmolog- 
ical data yields limited constraints on the scenarios in 
question, the data on the Cosmic Microwave Background 
(CMB) anisotropics expected in the next several years 
from the Planck satellite mission (scheduled to launch in 
2008) will lead to qualitative improvements. 

There are many models which generate strong neutrino 
interactions. A sterile neutrino added to the standard 
model may have its mass generated by a scalar field in 
analogy with the Higgs mechanism (though the mass- 
generating scalar is not, by necessity, the Higgs boson 
since the sterile neutrino carries no standard model in- 
teractions). That scalar will mediate neutrino interac- 
tions, and the size of the coupling is often large enough 
to make them measurable in the CMB. Cases in point 
are Mass- Varying Neutrinos (MaVaNs) [l^, [111 , a theory 
of neutrino dark energy, or a model of "late" neutrino 
masses [2^, [l^]. These theories feature neutrino-scalar 
interactions of a size which, as we will see, may be large 
enough to be measurable by the Planck satellite. We 
will see that Planck alone will be able to observe or rule 
out scenario with a single interacting and two standard 
neutrinos at the 4.2 u level. 

In many models, additional scalars may also become 
thermally populated. It is possible to choose the pa- 
rameters such that the scalars recouple after the BBN 
era [2^, [l^ . Hence an independent CMB constraint be- 
comes of great interest. Though the present constraints 
on Nl^^ are not very strong in comparison to the BBN 
constraints, we will also find that Planck will be able to 
observe a single extra neutrino at approximately 2.5 a 
level. 

It has been shown that in the case of interactions me- 
diated by heavy particles cosmology does not yield ef- 
fective new bounds [1^. Why is the situation different 
when the new particles are light? The basic idea can be 
easily stated. Coupling neutrinos to a light scalar field 
can populate additional degrees of freedom in the early 
universe, changing the expansion rate as well as the evo- 
lution of inhomogeneities. Moreover, it can drastically 
increase the neutrino-neutrino interaction cross section, 
much more than what is possible with heavy new par- 
ticles, turning the neutrinos in the early universe into 



a self-coupled fluid. The neutrino density perturbations 
would evolve differently from the standard case - coupled 
neutrinos would undergo acoustic oscillations rather than 
stream freely in the CMB epoch (T ~ 1 — 10 eV) - and 
affect the baryon-photon fluid and dark matter through 
their gravity during the epoch of radiation domination. 

Our cosmological analysis is done in a rather model- 
independent way. For most of the parameter space, 
one can simply assume that in the CMB epoch there 
is a certain number of free streaming relativistic species, 
Nps (the standard neutrinos fall into this category), and 
a certain number of relativistic species forming a self- 
coupled fluid, Ncoupied- This approach is certainly not 
new [2^, [s^l- We will denote the total number of ef- 
fective neutrinos (which includes any thermalized radia- 
tion, whether free-streaming or coupled) as N^^^. The 
results then have a broader applicability. For example, if 
one sets Ncoupied = 0, ouc obtains a bound on the total 
number of the standard neutrinos. The constraints on 
the latter from the current data, as well as forecasts for 
Planck, have been investigated in numerous recent stud- 
ies 0, [si m Hi, [33, 34, 3|, with somewhat differing 
results. We will weigh in on this controversy. 

The cosmology of neutrinos coupled to a light scalar is 
a very old topic. Over the years, different aspects of it 
were discussed in [H, 113, S, HI] , [H UlLFour modern 
data analyses are those by Hannestad [40|, by Trotta 
and Melchiorri [4l|, by Bell, Pierpaoli, and Sigurdson 
fsS'l, and by Cirelli and Strumia ^^V\. These papers only 
analyze the existing data and do not forecast the reach of 
Planck. Moreover, once again, the results these papers 
find with the available data differ. Wc will reexamine 
this situation using up-to-date cosmological data. 

In the next Section we examine in more detail the types 
of neutrino models which generate strong neutrino in- 
teractions or extra relativistic degrees of freedom. We 
then turn to the CMB phenomenology of constraining 
these models. We carry out an analysis using all the 
current cosmological data, and determine the precision 
with which Planck will be able to constrain neutrino in- 
teractions. We then apply these constraints to Majorana 
neutrino interactions and consider the implications for 
theories of neutrino dark energy (MaVaNs). The details 
of our analysis are contained in Sects. III-V; the Reader 
with interests in cosmology may wish to focus on these 
Sections. Conversely, the Reader who is interested only 
in the implications for neutrino models can proceed di- 
rectly from Sect. II to Sect. VI. 



II. MODELS OF NEUTRINO-SCALAR 
INTERACTIONS 

As mentioned in the Introduction, the possibility of 
coupling the neutrinos to a novel light scalar field has 
been entertained in the literature since the early 1980's. 
In this Section, we will briefly review some of the cosmo- 
logically relevant features of the many models that were 
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constructed. For more details, the Reader should con- 
sult, e.g., BSHIl,!!!!!!,!!!,!!!] and many other papers 
dedicated to the subject. 

From the point of view of cosmological constraints 
based on the CMB, we are interested in the physics at 1 
eV energies. At these low energies the relevant properties 
of the full models are captured by effective low-energy La- 
grangian terms. We may thus try to learn about some 
of the common features of the models by building this 
interaction "bottom-up" . 

The starting step may be to seek an interaction of the 
Yukawa type, g4>vv. To have the correct gauge struc- 
ture, the neutrino field must be promoted to the lepton 
doublet. At the level of dimension 4 operators one can 
then write L'^t • <pL. This means the new scalar field 
</> is a triplet of SU{2) and couples to the Z boson. If 
one further writes a symmetry breaking potential for 
(to obtain a massless Goldstone boson - a Majoron), one 
obtains the classical model of Gelmini and Roncadelli 
[§| . This model has been ruled out by the LEP measure- 
ments of the Z boson width: the triplet Majoron would 
contribute an equivalent of two extra neutrino species to 
the width. This argument extends to other models in 
which the new light scalar is not a singlet of SU{2); for 
instance, the Higgs doublets [HI, would contribute a 
half of an extra neutrino species to the Z width. As a 
consequence, we must assume that the light scalar field 
is a Standard Model singlet. 

The simplest renormalizable model that generates 
an effective low-energy neutrino-scalar coupling involves 
adding a right-handed sterile neutrino N to the theory 
and coupling a scalar </) to it. To this end, consider the 
standard see-saw mechanism. The mixing of N with the 
Standard Model neutrinos only enters through a Dirac 
neutrino mass term, so that the neutrino Lagrangian con- 
tains 

5C = yLHN + uinNN. (1) 

The first term is of the same form as the lepton Dirac 
masses in the Standard Model, while the second one is a 
Majorana mass term for the right-handed (sterile) neu- 
trino. When the sterile neutrino is integrated out, as- 
suming 771 jv S> = yv, only a light Majorana neutrino 
remains, with mass y^v^ jrnjq . 

We may now promote the sterile neutrino mass, m^r, 
to a dynamical (complex) field (/), in analogy to the Stan- 
dard Model Higgs mechanism, where the Dirac mass is 
generated by the dynamical Higgs [111 : 

8H = yLHN + X(t)NN. (2) 

When (j) develops a VEV, {(j>) — /, the sterile neutrino 
gets a mass tun = Xf and can be integrated out of the 
theory. The hght Goldstone mode G (</> = {f + p)e-''^/ 
is coupled to the light neutrinos via an effective interac- 
tion of the form 

^e^^lf{L{H)){L{H)). (3) 



Clearly, at the energies of 0(1) eV the Higgs field in 
Eq. ^ has no excitations, only the vacuum expectation 
values (VEVs), (H) = v. 

An alternative form of the interaction can be obtained 
by absorbing e^^^^ into the phase of L. The interaction 
then appears from the kinetic term for L, 

ydf^GLa^L + h.c. (4) 

This form makes manifest the derivative nature of the 
coupling of G. 

We note that while the construction reviewed here is 
the simplest, it is certainly not unique. In general, the 
scalar field need not couple to only the sterile neutrino 
and could be coupled via both the Dirac and Majorana 
terms in Eq. ([!]), e.g., by promoting y 0/A, tojv 
$VA [H. 

Expanding Eq. ([3]) to the first order in G, we find a 
Majorana mass for L, rrii, = y'^v^ /{\f) (where v = {H)) 
and an effective Yukawa (pseudoscalar) interaction of G 
with L, with the coupling strength 

_ y'^v'^ _ Am|, _ my 

The second equality has a very simple physical interpre- 
tation: rnjj/mj^ gives the amount of admixture of the 
sterile neutrino into the light mass eigenstate and A is 
the coupling of g to the sterile component. 

In order to generate a sizable coupling g between the 
neutrino and scalar, there must not be too large a hier- 
archy between rrii, and the scale of symmetry breaking /. 
If the sterile neutrino mass is at GUT scale, the coupling 
of G is of the size g ^ 10~^^, assuming ~ 0.1 eV 
and A ~ 1, which is too small to give rise to significant 
effects during the CMB epoch. On the other hand, if the 
scale of new physics, A, is TeV scale, couplings of the size 
g ~ lO"'^'^ may be generated. Though this is a very small 
coupling, it is large enough to result in observable effects 
through decay-inverse decay of the neutrino to scalars, 
which tightly couples the scalars and neutrinos, remov- 
ing neutrino free-streaming. We return in detail to the 
constraints on the coupling g in Sect. I VII 

There are many theories where a coupling g > 10"^'^ 
is generated. Theories of neutrino dark energy introduce 
interactions in the effective theory of the form Eq.[21 often 
with a sizable couplin g g. In the models of dark energy 
originally discussed in [43| | and revived in [i^ as neutrino 
dark energy, the scalar generating the dark energy is a 
pNGB. If the scale of the symmetry breaking / of the 
scalar U(l) symmetry is sufficiently low (TeV scale or 
lower) , a cosmologically interesting coupling between the 
scalar and neutrino may be generated. Since the scalar is 
a pNGB in these models, it is protected from large radia- 
tive corrections, and hence remains light-an advantage 
in any theory of dark energy. In a similar model of [4q . 
/ ~ Mpi, but mpf and ttid are small, through extremely 
small Yukawa couplings y and A to the neutrino mass 
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generating fields H and <j). Since A is very small, g is also 
very small and there are no observable consequences in 
the CMB through the neutrino-scalar coupling. Because 
the sterile neutrino is light, however, and the mixing be- 
tween the active and sterile neutrino is substantial, the 
sterile neutrino may become thermalized, and there may 
be signals in the CMB through increased Nps- 

Another theory of dark energy may also generate in- 
teresting signals in the CMB. In contrast to these mod- 
els where the scalar field is typically associated with 
a broken symmetry at a comparatively high scale /, 
MaVaNs, as introduced in [23|, place the entire sector 
around a meV-this scale is chosen according to the mea- 
sured neutrino mass splittings and the dark energy scale 
5m1 ^ 10~^ eV, A^b ~ 10~^-^ eV. As a consequence, 
all mass scales in the model, including the Dirac mass 
and the scalar mass, lie in the sub-eV range, and the 
coupling A is typically not too small; the small hierarchy 
between the Dirac and sterile neutrinos implies a sizable 
coupling g in many cases. 

There are many other instances where the hierarchy 
between the Dirac neutrino and sterile neutrino is much 
smaller due to the introduction of lighter sterile neutri- 
nos. Much lighter sterile neutrinos have been considered 
in a wide variety of contexts, most notably perhaps in 
connection with the LSND measurement, where the pres- 
ence of a sterile neutrino with mass around 1 eV has been 
invoked to explain the appearance of muon neutrinos [i^ . 
Other models feature a keV mass sterile neutrino as the 
dark matter [13, l48l. 49| , sometimes with accompanying 
keV mass scalars [50 1, and weak scale neutrinos associ- 
ated with SO(IO) GUTS, where the addition of the neu- 
trino is necessary for anomaly cancellation [sij. Various 
low energy see-saws have also been considered, as in [26| . 
where the scale / is in the 50 MeV to 500 GeV range. 

We have seen that there is a broad class of models 
which generate exotic neutrino-scalar interactions. The 
addition of these interactions, depending on the choices of 
X, mo and m^r, may be good candidates for observation 
in the cosmic microwave background. First, additional 
scalars may become thermalized and increase the effec- 
tive number of neutrino species. Second, these scalars 
mediate additional neutrino interactions, through scalar 
mediated neutrino scattering, which could remove neu- 
trino free-streaming at CMB temperatures. We now turn 
our attention to studying the impact of non-standard be- 
havior in the dark radiation sector on CMB anisotropies. 



to consider both the effects of the neutrino background 
density and that of the neutrino density perturbations. 
Nl;f^ dictates the expansion rate of the universe and, to- 
gether with matter density fim^^, controls the redshift of 
matter-radiation equality, z^q. The latter effect is more 
subtle. Assuming adiabatic initial conditions, the neu- 
trino and photon inhomogeneities are of comparable size 
to begin, so the presence of neutrinos modifies the evo- 
lution of the photon perturbations. When a perturba- 
tion of a given size enters the horizon, the gravity of 
the neutrino perturbation is comparable to the gravity 
of the photon perturbation. The subsequent evolution 
of the two, however, are different. The photon-baryon 
plasma oscillates like a compressible fluid; the standard 
neutrinos, on the other hand, stream freely, quickly eras- 
ing their density fluctuations. Gravitational coupling be- 
tween the two means that the evolution of fluctuations in 
the photons could be affected by neutrino free-streaming. 
In fact, Stewart 52j noted back in 1972 that if this effect 
were large, it would jeopardize structure formation. 

Shortly thereafter, Peebles [5^ numerically solved the 
problem of the coupled evolution of the neutrino and pho- 
ton perturbations. He showed that the neutrino inho- 
mogeneities in the standard case do indeed decay shortly 
after entering the horizon; in the process, they damp pho- 
ton inhomogeneities, although at a significantly smaller 
level ('^12%) than anticipated by Stewart. Much more 
recently, the problem was addressed by more accurate 
numerical computations and analytically [3^ . [s^ l and for 
the CMB power spectrum the amount of suppression (rel- 
ative to the case where all the radiation is strongly cou- 
pled) was found to be 



SCi/Ct ~ -0-53 PFS/Prad, 



(6) 



where pps is the energy density in freely streaming ra- 
diation and prad is the total radiation density (free- 
streaming plus strongly coupled, including photons). 

Coupling the oscillating photons to the damped neu- 
trinos through the gravitational potential also changes 
the phase of the photon oscillations, due to a change in 
the speed at which information propagates in the fluid 
from the presence of freely streaming neutrinos, com- 
pared to the hypothetical case of no free-streaming neu- 
trinos. This effect was clearly established in f33|. The 
resulting shift of the CMB peaks is 
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III. TIGHTLY COUPLED NEUTRINOS: 
MODIFIED EVOLUTION EQUATIONS. 

We summarize the relevant physical effects of dark ra- 
diation (i.e. neutrinos) on the CMB. The energy density 
in relativistic neutrinos [qJ] is a fraction ~ N^^ ^ / {N^^ ^ + 
4.4) of the total radiation (freely-streaming and cou- 
pled neutrinos plus photons). Thus during the radia- 
tion epoch the neutrino gravity is important. One needs 



Remarkably, both the amplitude suppression and 
phase shift are clearly present in the Peebles' solution 
[95j . These effects are at the core of the physics behind 
the sensitivity of CMB to neutrino free-streaming. The 
amount of damping and the phase shift change if either 
additional freely streaming relativistic ( "neutrino- like" ) 
species are added or the neutrinos become self-coupled. 
In the latter case, both effects are removed: the neu- 
trino fluid oscillates similarly to the photon fluid (with- 
out baryon loading). 
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Notice that the effect is not uniform for all CMB mul- 
tipoles. The above argument was made for modes enter- 
ing the horizon in the radiation era. For modes entering 
the horizon in the matter era, the neutrino perturbations 
have no effect. Thus, the damping is operational only on 
small scales, I > 200. 

In our analysis, we consider interacting neutrinos in 
the tightly coupled limit. By this we mean that the neu- 
trinos can be approximated by a fluid for the entire range 
of relevant scales, including the scales corresponding to 
the CMB multipoles of / ^ 2000 and those measured by 
the large scale structure (LSS) surveys SDSS and 2dF. 
The neutrino analogue of the Silk damping scale is thus 
assumed to be well below ^ 0(10) comoving Mpc. The 
implications of this assumption are further discussed in 
Sect.lVll 

With the above assumptions, the Boltzmann equations 
for the coupled neutrinos are very simple, as discussed in 
[55| : the standard multipole expansion for the neutrino 
perturbations (see [sl]) is truncated at the level of density 
and velocity perturbations. The quadrupole (shear) and 
higher order moments of the perturbations are set to zero. 
The analogues of Eqs. (49) or (50) in [5^ are: 

• Synchronous gauge 
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(J^ = 0. 
• Conformal Newtonian gauge 



(8) 
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Here all the conventions arc those of [56|. The quan- 
tity Si, = Spiy/pi, is the neutrino density perturbation; 
01, = ik^ ST'^j / {p + P), with the "i^" index assumed on the 
right hand side, is the neutrino velocity perturbation; a 
is the shear (see Eq. (22) of [1^). The quantity h is one 
of the two scalar perturbations in the synchronous gauge 
(the one corresponding to the trace of the scalar metric 
perturbation). and ip are the scalar metric perturba- 
tions in the Conformal Newtonian gauge. They coincide, 
up^to a sign, with the gauge invariant Bardeen variables 

The same method of truncating the multipole expan- 
sion is utilized in the earlier analyses. Our equations 
agree with those of [s^, H^l (which are stated in the syn- 
chronous gauge) and those of (sTj] (which are stated in the 
Conformal Newtonian gauge). We note that the CAMB 
[sst code we use (also CMBFast [11]) employs the syn- 
chronous gauge. A slightly more general parameteriza- 
tion in terms of the "viscosity parameter", c^^^ [6(| is 
followed in lilll. 
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FIG. 1: Effect of neutrino free-streaming on the CMB multi- 
pole spectrum. The thickest curve is the spectrum with the 
best fit parameters from WMAP3; the other curves (from 
bottom to top) correspond to 1, 2, and 3 strongly coupled 
neutrinos (keeping the total number of neutrinos fixed at 3). 



The effect of making the neutrinos coupled, for fixed 
cosmological parameters, is illustrated in Fig. [1] The 
thick curve refers to the standard case of three freely 
streaming neutrinos, while the other curves illustrate the 
effect of coupling 1, 2, and 3 neutrinos (in order of de- 
viation from the thick curve) . The changes of the tem- 
perature (TT), polarization (EE), and the cross correla- 
tion between them (TE) are shown. The Figure clearly 
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exhibits both the amplitude suppression and the phase 
shift effects at I > 200. 

The difference in the anisotropics on scales / > 200 
is quite large (on average about 25-30% for three cou- 
pled neutrinos), well outside the errors of the current 
WMAP data. However, this does not mean than the 
coupled neutrino scenario is already excluded. Indeed, 
it must be kept in mind that there are many cosmolog- 
ical ("nuisance") parameters that can be adjusted, such 
as the baryon, dark matter and dark energy densities, 
the spectrum of primordial fluctuations, and others. By 
adjusting these parameters, it may be possible to undo 
most of the effect of the neutrino self-coupling. 

This issue of "degeneracies" between different param- 
eters is of course well known in cosmology. A sim- 
ple illustration of it is given in Fig. [2] where we show 
changing the CMB power spectra as one changes the to- 
tal number of freely streaming neutrinos. While sim- 
ply changing Nps to 7 produces a large shift in the 
position of the peaks (because of the faster expansion 
in the radiation era as discussed above), the effect can 
be compensated by changing other parameters such that 
the redshift of equality is preserved 1 + Zeq — 

4.05x10^^2™ /iV(l+0.6905iV^^//3.04) is fixed while vary- 
ing N^-f^ by scaling to compensate for the increase in 
N^-f^ , while also fixing the physical baryon density flbh^ 
and flm- Indeed, the physical quantities that are mea- 
sured in CMB are dimensionless quantities (angles on 
the sky), hence they depend on the ratios of the physical 
densities, etc. See, e.g., [6l]] for further discussion. 

Not all effects follow this simple argument. For exam- 
ple, the Silk damping does not, as it involves a physical 
dimensionful constant, the Thompson cross section. The 
faster expansion of the universe with more neutrinos im- 
plies more Silk damping in the high multipoles of the 
CMB j36j. This can be partially compensated by adjust- 
ing the Helium fraction, as discussed in [s^l . It should be 
kept in mind that this mechanism is limited by a variety 
of astrophysical considerations. We will return to this 
topic in Sect. IV Bl 

The real challenge is to establish the size of the resid- 
ual differences of the CMB predictions in the two sce- 
narios, after appropriately adjusting the "nuisance" pa- 
rameters, in comparison with the resolution of the exper- 
iments. These residual differences turn out to be much 
smaller than the differences seen in Fig. [T] This fact ren- 
ders difficult writing down a simple estimate for the pre- 
dicted sensitivity of the present and future experiments 
using order-of-magnitude arguments and necessitates a 
detailed scan of the multidimensional parameter space. 

In Sect. IV Bl we show how well the effects of amplitude 
suppression ^ and phase shift ([7]) can be compensated 
by adjusting the cosmological parameters and how big 
the residual differences are. We will also see which CMB 
multipoles are essential for testing the neutrino sector 
and how robust our predictions for Planck are. A com- 
plete analysis of this type has not been done before. 

We now present the results of our numerical studies. 
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FIG. 2: Effect of extra neutrinos on the CMB multipole spec- 
trum. The central black curve is the spectrum with the best fit 
parameters from WMAP3, the top green curve the spectrum 
with 7 freely streaming neutrinos, and the lowermost magenta 
curve results when the total number of freely streaming neu- 
trinos is 7, but Zeq is kept fixed by varying . 
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IV. NUMERICAL RESULTS 
A. Sensitivity of the current data 

1. Literature overview 

We begin by summarizing the bounds on the numbers 
of freely streaming and self-coupled relativistic species 
derived in the literature. In the case of extra relativis- 
tic species, a number of studies have been undertaken in 
the last several years. Crotty, Lesgourgues, and Pastor 
found the effective number of neutrinos to lie in the inter- 
val [1.4, 6.8] at 95% confidence level (C.L.), by combining 
WMAP year one data |i62|] with the LSS data from 2dF 
[11] and a prior on the Hubble constant. Similar results 
were found by Pierpaoli [63l and Hannestad [65^ . Bell, 
Pierpaoli, and Sigurdson [30| have shown that WMAP 
year one data [62| , Lyman-a forest [i^ , and Sloan Digital 
Sky Survey (SDSS) data [gt'I taken together constrain the 
effective number of neutrinos to be 4.9I2 3 ^5% C.L. 
Seljak, Slosar and McDonald y| have also computed this 
number with more recent data - including the WMAP 
three year data (WMAP3) ^6^, i63|, Lyman-a forest data 
[to! ]. SDSS measurements of the baryon acoustic oscil- 
lations [7l[ (see for a complete list). They find a 
similar result, Nps = 5.3li;7 at 95% C.L. The WMAP 
team finds year-3 data combined with large-scale struc- 
ture and supernovae prefer 3.3±1.7 neutrinos 72]. Cirelli 
and Strumia [3l[ find Nps = 5 ± 2 (2tT errors). In con- 
trast, Hannestad and Raffelt [3^ . constraining neutrino 
mass and relativistic energy density together, determine 
2.7 < Nfs < 4.6 at 95% C.L. using WMAP3 and large- 
scale structure data from 2dF and SDSS. 

In the case of self-interacting neutrinos, all fits are 
consistent with no interactions, but with (very) differ- 
ent confidence levels. Hannestad [13] concludes that in 
the model with three interacting neutrinos "it is impossi- 
ble to simultaneously fit CMB and LSS data" . From his 
Fig. 8 we read off a 3a exclusion using only WMAPl and 
6a exclusion when this data is combined with the SDSS 
LSS data and the HST Hubble constant measurement (in 
the limit of zero neutrino mass). In contrast, Trotta and 
Melchiorri [4l| find that the interacting neutrino scenario 
is disfavored at only 2Aa. Similarly, Bell, Pierpaoli, and 
Sigurdson find Ncoupied < 3.0 at 95% C.L. CirelH and 
Strumia find the number of additional interacting neu- 
trinos constrained to be less than 1.3 at la with more 
recent data (including WMAP3). 

In short, there is quite a bit of variation between the 
published results. Some of this variation could be at- 
tributed to different data used in the calculations (for ex- 
ample, WMAP 1-year vs. WMAP 3-year data releases, 
whether Lyman-a was included in the analysis or not, 
etc), but some clearly must be due to the differences in 
the analyses themselves. We therefore consider it well- 
motivated to repeat the calculations, using the most re- 
cent data available to us. In Sect. IIV A 21 we consider the 
sensitivity of WMAP3 alone, while in Sects. IIV A 31 and 



IIV A 41 we add other cosmological data. 

As far as forecasting for Planck, published results on 
Nps likewise differ. In addition, no analysis on Ncoupied 
has been performed. Ref. [33] finds that Planck will be 
able to constrain the total number of freely streaming 
neutrino species to ANps — 0.24 {la error) using tem- 
perature information only and a sky coverage /sky = 0.5. 
In contrast, Ref. 34] finds using Planck temperature in- 
formation only ANps = 0.6 even with a more optimistic 
sky coverage /sky = 0.8. Both studies employ the same 
technique (the Fisher matrix analysis). Finally, Ref. [35] 
investigated the bounds on Nps using both the Fisher 
analysis and Markov Chain Monte Carlo (MCMC). The 
results are quite intriguing: while the Fisher analysis 
yields ANps = 0.26(0.27) with (without) lensing, the 
MCMC method yields very different results depending 
upon whether lensing is assumed - with lensing the re- 
sult is consistent with the Fisher analysis, while without 
lensing ANps = 0.46. The forecast for Planck's sensitiv- 
ity to the number of self-coupled neutrinos to the best of 
our knowledge has not been done. We perform a com- 
bined analysis of Planck's sensitivity to freely streaming 
and self-coupled neutrino species in Sect. IIV B II We also 
investigate, in Sect. IIV B 21 a potential impact of combin- 
ing Planck with other cosmological data. 



2. WMAP3 alone 

As a first calculation, we investigate how sensitive 
WMAP is by itself to neutrino self-coupling. We do this 
by fitting WMAP3 [H,!!^, varying a set of cosmological 
parameters (to be specified shortly) under four different 
assumptions: (i) three standard free-streaming neutri- 
nos; (ii) one neutrino coupled, two free-streaming; (iii) 
two neutrinos coupled, one free-streaming; (iv) all three 
neutrinos coupled. We then compare the goodness of 
fit at the corresponding best-fit points. We also explore 
the sensitivity of WMAP to the total number of (stan- 
dard) neutrinos, in a similar way. For this, we consider 
two more scenarios: one with the total neutrino number 
equal to five and another with a single neutrino flavor. 

The fitting is done using the MCMC code COSMOMC 
[Zl, with the CAMB code [H] modified by us to in- 
clude both freely streaming and tightly self-coupled neu- 
trinos. The MCMC method is by now a standard tool 
in cosmology, used for both data analysis and forecast- 
ing. It is employed by the WMAP [tI, [tJ and Planck 
[75| teams, as well as many other groups (in particular, 
among the papers reviewed in Sect. IIV A li [2„,30., ,35, ,4i| 
use MCMC). With the MCMC method, the likelihood 
function is mapped out in a multidimensional region of 
parameters around its maximum. As a result, one gets 
the location of the best-fit point, with the correspond- 
ing likelihood characterizing the goodness of fit, as well 
as the allowed ranges of the parameters. One need not 
make a priori assumptions on the functional form of the 
likelihood function, although choosing the parameteriza- 
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TABLE I: The best fit values of —log(likelihood) (second column) and the corresponding shift of the best fit parameters as 
a function of the number of coupled neutrinos. The last five columns list derived parameters. The initial spectrum is defined 
with a pivot point kpiv ~ 0.05 Mpc~^. The fit is to the WMAP3 dataset only. 



tion in such a way that the posterior distributions are 
approximately Gaussian and there are no strong correla- 
tions saves computer time Q [97j . 

In this Subsection, we choose the following six cos- 
mological parameters: the physical densities of baryons, 
flbh^, and dark matter, flch^, the ratio of the approxi- 
mate sound horizon to the angular diameter distance 9, 
the optical depth to the last scattering surface r, and 
the primordial spectrum of the scalar curvature pertur- 
bations, characterized by the spectral index and the 
power As on a preset ("pivot") scale (taken in this cal- 
culation to be fcpTO = 0.05 Mpc~^). We assume that the 
Universe is flat, there is no running of the spectral in- 
dex, the Helium fraction is fixed at Yho — 0.24 and the 
dark energy is a cosmological constant {w = —1). As we 
will see, in the case of WMAP3, the six parameters we 
vary contain the necessary degeneracies with the neutrino 
coupling. 

The results are tabulated in Table HI The first row cor- 
responds to the standard scenario, Nps = 3, Ncoupled — 
0. The next three rows correspond to coupling one, two, 
or three of the neutrino species, while keeping the total 
number of neutrinos at three. Finally, in the last two 
rows, we vary the total number of neutrino species, as- 
suming they are all freely streaming. 

The second column shows the difference between the 
of the best fit in a given scenario and the correspond- 
ing quantity in the standard scenario (first row). The 
third column shows the corresponding confidence level 
(C.L.) [ol]. The next six columns in the Table show 
the values of the best-fit parameters for each scenario. 
The last five columns show the values of derived param- 
eters: the cosmological constant Ha, the age of the Uni- 
verse, the matter fraction Qrm the red-shift of ioniza- 
tion Zre and finally the Hubble constant, Hq. To com- 
pensate for the effect of neutrino coupling, Hq increases, 
while Qm and the spectral index decrease. To compen- 
sate for additional neutrino species, all three quantities 
increase. In the latter case, it is easy to check that the 
best-fit parameters shift in such a way that the redshift of 
matter-radiation equality, Zeq, is preserved, as discussed 
in Sect. IIIII As also increases with Nps to compensate 
for the reduction of perturbations on all scales. 

From the second column, we see that the effects of the 
neutrino coupling as well as the variation of the total 
number of neutrinos are nearly perfectly compensated. 



The quality of the fits is virtually the same in each of the 
six cases. Therefore, WMAP by itself cannot distinguish 
between these scenarios. This conclusion is consistent 
with the findings of the analysis by the WMAP collabo- 
ration, which considers the sensitivity of the experiment 
to the number of freely streaming neutrino species [t^ . 
It differs from the findings of _J£(| where a Scr exclusion 
of the {Nps = 0, Ncoupled = 3) scenario was claimed. 
Clearly, our conclusion relies on the ability of the code 
to find the set of parameters that give the most complete 
compensation of the effects of modifying the neutrino sec- 
tor. We see that the MCMC method accomplishes this 
task well. 



3. WMAP3 plus large scale structure (2dF, SDSS), HST, 
SN la. 

We now explore if the addition of the other presently 
available cosmological data can lift the degeneracy of 
WMAP3. In this analysis, we include the data on the 
large scale structure (LSS) from the 2dF [t^ and the 
SDSS survey - both the Main l&f] and the recently re- 
leased Large Red Galaxy (LRG) [77] data samples. We 
also use the Type la supernova data from the SNLS col- 
laboration [78], as well as the Hubble constant measure- 
ments from the Hubble Space Telescope (HST) [t^, in 
the form of a Gaussian prior h — 0.72 ± 0.8. 

We vary eight parameters, six as described in the pre- 
vious Subsection, plus two more, the helium fraction, 
Ihc and the running of the spectral index, Ug. These 
two parameters are added in anticipation of the analy- 
sis for Planck, where their roles are important. With 
the present day data, as we will see, they do not make 
much of a difference. For the helium fraction, we use a 
gaussian prior Yhc — 0.24 ± 0.009 motivated by obser- 
vational bounds [80|. We also impose a rather generous 
hard prior, 0.15 < Yne < 0.3. The values greater than 
0.3 would be problematic for the solar model (since the 
input value into the solar model cannot be less than the 
primordial value) [81, 83]. 

We begin by performing a combined fit to all the data. 
Results for several specific scenarios are shown in Ta- 
ble |TT1 We find that the de generacies of the WMAP3-only 
analysis are broken by the additional data. In particu- 
lar, relative to the standard case of three freely streaming 
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TABLE II: The same as Table HI but including WMAP3, SDSS Main and LRG data samples, 2dF, HST, and SN la datasets. 
The last column is the derived best fit value of as. 
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FIG. 3: The shift of the best-fit parameters as a result of neutrino coupling. The dashed-dotted curves correspond to the 
scenario of coupled neutrinos, while the solid ones refer to the standard (free streaming) neutrinos. In both cases, the fits are 
to the datasets from WMAP3, SDSS, 2dF, HST, and SN la. 



neutrinos, the case with all three neutrinos coupled is dis- 
favored at 3.5(7, while two and one coupled neutrinos are 
disfavored at 2.5a and 1.5a, respectively. Assuming no 
self-coupled neutrinos, the scenario with a single freely 
streaming neutrino is disfavored at 3.1a, while the sce- 
nario with five freely streaming neutrinos gives a fit which 
is just as good as the standard one. 

The corresponding best-fit parameters are tabulated in 



Table |TT1 One can see that the general trends are sim- 
ilar to what was observed with WMAP3 only: as more 
neutrinos are coupled, the fit prefers larger values of the 
Hubble constant and smaller values of the spectral in- 
dex Ug. The shifts are graphically illustrated in Fig. [31 
The standard free streaming neutrinos are shown in solid 
curves, while the scenario of three self-coupled neutrinos 
is shown with the dashed-dotted curves. 
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FIG. 4: The combined reach of all available data (excluding Lyman-a) on the number of free-streaming neutrinos, Nfs, and 
the number of coupled neutrinos, N^oupUd.: WMAP3 + SDSS (LRG) + SDSS (Main) + 2dF + HST + SNIa. The sohd contours 
indicate 1 and 2a C.L. 



From the point of view of particle physics one 
may be interested to know how the number of al- 
lowed self-coupled neutrino species varies with the 
total number of neutrinos. As mentioned earlier, 
one can easily imagine models in which additional 
neutrino-like degrees of freedom are populated after 
the time of the BBN but before CMB (see, e.g., 
[27II). In Fig. [4] we present the allowed region in the 
{N F s 1 N coupled) plane. The plot was obtained by run- 
ning MCMC in the ten-dimensional space of parameters 
{Nfs, N coupled, ^bh'^,^eh'^,0, r, Yhc 

, fig , ^run ; 

log[10i%]), 

assuming a ACDM universe, and marginalizing over all 
but the first two parameters. From the Figure, we see 
that scenarios with no freely streaming neutrino are 
clearly disfavored, whether the number of self-coupled 
neutrinos is zero or three. At the same time, the con- 
straint on Ncoupied is relaxed if in addition to the coupled 
neutrinos there are also ^3 — 5 freely streaming ones. 
In fact, even with standard neutrinos the fit actually 
prefers the total number of freely streaming neutrinos 
to be greater than three: the best fit is achieved for 
Nfs = 3.7, Ncoupied = 0. This curious result warrants 
further investigation. 

It turns out that the piece of data responsible for fa- 
voring large Nps is the Main data sample from SDSS. 
With it removed (i.e., with the matter power given by 
the LRG SDSS dataset and 2dF dataset), the fit changes 
dramatically, as shown in Fig. [51 The new fit has a slight 
preference for values of Nfs ^ess than three: the best fit 



lies at Nfs = 2.5, Ncoupied = 0.1. 

Conversely, if we remove the LRG data from the fit, 
the best fit point moves all the way to Nfs — 5.2, 
Ncoupied — 0.1, as can be seen in Fig. [SI Marginahz- 
ing over all other parameters including Ncoupied yields 
Nfs = 4li:3. Thus, it may be too early to conclude that 
the SDSS LRG and Main samples are consistent with 
each other. At least as implemented in COSMOMC, they 
pull the best-fit value of Nfs in different directions, with 
the LRG sample favoring the standard values. 

Our results show that if one chooses to rely on the 
LRG sample, one finds much less sensitivity to neutrino 
self-coupling: the point Nps — 0, Ncoupied — 3 lies inside 
the 2 a contour, while Nps = 1, Ncoupied = 2 lies on 
the 1 a contour. Thus, at present the bounds on coupled 
neutrinos from the global fit should perhaps be taken 
with caution. 



4. Adding the Lyman-a data. 

We also repeat the analysis including the Lyman-a 
dataset in the fit. This dataset has been somewhat con- 
troversial [9^ . Our main reason for doing this calculation 
is to get an idea about the additional sensitivity that can 
be gained from this dataset, and also to further check 
consistency with the published analyses where this data 
is used. 

We find that the overall sensitivity to coupled neu- 




FIG. 6: Sensitivity of WMAP3 + LSS from SDSS and 2dF + HST + SNIa. Compared to the previous figure, LRG data has 
been replaced with SDSS Main sample. SDSS main prefers larger Nfs compared to LSS from LRG data sample. 
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trinos is somewhat increased with the addition of this 
piece of data. Relative to the standard case of three 
free-streaming neutrinos, the case with a single coupled 
neutrino is disfavored at l.Scr, the case with two coupled 
neutrinos is disfavored at 3.1a and, finally, the case with 
three coupled neutrinos is disfavored at 4.1(t. In other 
words, the addition of the Lyman-a data does bring fur- 
ther improvements in sensitivity, though not very large 
ones. The Lyman-a dataset does favor large values of 
A^F?, just like the SDSS Main sample, as also observed 
in [20] , • The comments at the end of the last Subsec- 
tion apply here as well. 



B. Sensitivity of Planck 

1. Planck only 

As already mentioned, the situation will improve 
markedly with the expected data from the Planck satel- 
lite. In this subsection, we describe the result of our 
MCMC analysis of Planck's sensitivity. 

We generate mock data for Planck in the all_l_exact 
data format of COSMOMC, using the best-fit point of 
WMAP3 as the "true" (seed) value and assuming flat 
Universe, w — —1 for the cosmological constant, and 
Nps = 3.04, Ncoupied — for the numbers of neu- 
trinos. The characteristics of the Planck detector (the 
beam size and the noise levels for temperature and po- 
larization measurements) are given in [7^]. The rele- 
vant measurements will be done in three frequency chan- 
nels of the High Frequency Instrument (HFI), 100, 143, 
and 217 GHz. In the literature [H, [IJ], the analy- 
ses have been done assuming an effective single channel 
for Planck. Ref. [33| in particular uses the sky cover- 
age of fsky = 0.80, the detector noise for temperature 
Wrp^^"^ = 40/zK', and for polarization Wp^^^ — 56/iK' and 
the beam size 9b — T (see Sect. IV Al for the definitions of 
these quantities). We have checked that these effective 

numbers (especially the value oiwp^^'^) are in reasonably 
good agr eement with these three-channel parameters of 
[Tq [lO0j |. The only substantial difference is th at we be- 
lieve it is more appropriate to take fsky = 0.65 |lOl| . 

Our analysis here parallels that of the last Subsection. 
We first fit the mock data in the several specific sce- 
narios considered early, namely, varying the numbers of 
self-coupled neutrinos keeping the total neutrino number 
at three, and then varying the number of standard neu- 
trinos assuming no self-coupled neutrinos. We compare 
the quality of the best fits in each case and observe how 
the cosmological parameters change to compensate for 
the effects of neutrino coupling. For our second analy- 
sis, we perform a global fit in the space of parameters 

log[10io^]) 

We then marginalize over the last 8 parameters to obtain 
the allowed region in the space of Nps, iV™"P''^'^. 

The results of the first set of calculations are tabulated 



in Table IIIII The trends in the parameter shifts are sim- 
ilar to the case of WMAP3 - both Ug and Ag decrease 
while the physical baryon density increases as more neu- 
trinos are coupled - but with Planck the shifts are much 
constrained. Moreover, the degeneracies are very cfh- 
ciently broken and the quality of the best fits is signifi- 
cantly poorer in the coupled cases relative to the standard 
case. The scenario of a single coupled neutrino {Nps = 2, 
Ncoupied = 1) is disfavored at A.2a relative to the stan- 
dard one. The scenario with {Nps — 5, Ncoupied = 0) is 
disfavored at 5.4ct, while the other scenarios in the Table 
are ruled at greater than 8a. Clearly, Planck's sensitivity 
will be dramatically better than that of the current data. 

The situation is graphically illustrated in Fig. [71 where 
the solid curves show the expected measurements at 
Planck and the dashed curves are those obtained with 
the combined current data (see Sect. IIV A 3|) - both un- 
der assumption of three standard freely streaming neutri- 
nos. The dotted curves show how the best-fit parameters 
shift if one instead fits the data under the assumption 
of a single coupled neutrino {Nps = 2, Ncoupied — !)■ 
Clearly, Planck by itself will have errors that are signif- 
icantly smaller than those of today's experiments com- 
bined. 

The results of the second calculation show that this 
level of sensitivity persists for any direction in the 
(N PS, Ncoupied) plane. The sensitivity contours in this 
plane are depicted in Fig. [51 We wish to stress two main 
results of this analysis: (i) the sensitivities of Planck to 
Ncoupied and Nps are similar, and (ii) for no direction 
in the plane {Ncoupied, Nps) is there is a degeneracy be- 
tween the two parameters. We will return again to the 
last point in Section [VAI 

Our prediction for Planck's sensitivity to the number 
of freely streaming neutrinos is asymmetric, ANi^ =io 3- 
The lower error roughly agrees with that of [SJ] , when the 
choice of higher fsky is taken into account, and (ssj ob- 
tained with the Fisher matrix method. It is in reasonable 
agreement with what [3^ finds with the MCMC method 
assuming no lensing of the CMB. We have investigated 
the effect of lensing on the sensitivity and we do not find 
the large effect of [3^1 • Instead we find the bounds re- 
main similar, slightly weaker for Nps, slightly stronger 
for Ncoupied- Notice that 34], considering Planck's sen- 
sitivity to Nps with the Fisher analysis, also finds the 
effect of lensing small. 



2. Planck plus other cosmological data 

Lastly, we briefly consider what might be gained by 
combining Planck with the other cosmological data. For 
that, we run a combined fit to Planck, LRG and Main 
samples of SDSS, 2dF, the Type la supernova data and 
the Hubble constant measurements from the HST. 

The result of the fit is shown in Fig. [51 The sensitivity 
to the number of freely streaming neutrinos is unaffected 
by the addition of the other data, while the sensitivity 
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TABLE III: The same as in Table m but now fitting to the mock Planck dataset. 




FIG. 7: The improvement of sensitivity of Planck (solid) over the current data (dashed), assuming the standard scenario of 
three freely streaming neutrinos. The dotted curve illustrates how the parameters measured at Planck would shift assuming 
a single neutrino species is coupled. The current dataset is comprised of the data from WMAP3, SDSS Main and LRG, 2dF, 
HST, and SN la, plus the observational prior on Yhc (see text for details). 



to coupled neutrinos improves by a factor of two. This 
can also be seen from the one-dimensional marginahzed 
plots in Fig. [lOl 

One should, of course, not overinterpret this result. It 
should be kept in mind that the mock data for Planck 
was generated from the best-fit point of WMAP3 and we 
did not explore the dependence of the fit on this choice. 
Also, by the time Planck's data is available, the other 



experiments could be updated. Nevertheless, the general 
lesson is that while most of the precision will come from 
Planck, the addition of the other data may lead to further 
tightening of the bounds. 




FIG. 9: Sensitivity of Planck plus other cosmological data. The constraints on Nps are not much changed as compared with 
Planck alone, but the figure shows a factor of two potential improvement in Ncoupied- 
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FIG. 10: Sensitivity of Planck plus other cosmological data (dashed-dotted) vs. Planck only (solid). 1-d marginalized proba- 
bilities. 



V. DISCUSSION OF THE RESULTS AND 
FURTHER NUMERICAL INVESTIGATIONS. 

The results of the previous Section indicate that Planck 
should constrain the numbers of coupled (and freely 
streaming) neutrino species much more accurately than 
what is possible today. In this Section, we make this 
statement more robust by investigating the following 
points: 

• The sensitivity is an indication that not all of the 
effects of neutrino coupling could be compensated 
and that the residuals should be within the sen- 
sitivity reach of Planck's instruments. We recall 
that in the last Section, we varied these param- 
eters: NfS, Ncoupled, Wb, t^c, 8, T, Us, Umn, Ag, 

Yhc- We have omitted other cosmological parame- 
ters such as the curvature ^Icurv , or the dark energy 
equation of state w. If varying these additional pa- 
rameters introduced additional ways of compensat- 
ing the effects of neutrino coupling, the bounds of 
the previous Section could be weakened. We need 
to show that this is not the case. 

• We also made assumptions about the performance 
of Planck's detector. We need to check the robust- 
ness of our results with respect to changing the 
characteristics of the detectors, such as the angular 
resolution and the detector noise levels. Put an- 
other way, we need to establish which multipoles I 
are crucial for this measurement. 

We will investigate these points using the Fisher matrix 
technique. 

In addition to these practical issues, we will also in- 
vestigate the issue of principle: what is the nature of 
Planck's sensitivity to neutrino coupling? We do that by 
examining which effects of coupling cannot be compen- 
sated by adjusting the cosmological parameters. 



A. The role of other cosmological parameters 



technique. The idea is very simple: suppose that the 
likelihood function around the best-fit point is approxi- 
mately gaussian. In this case, one can fix the parameters 
of the gaussian (in the case of an N dimensional param- 
eter space, an N X N symmetric matrix) by evaluating 
the likelihood function at 0(20) points, rather than at 
10^ points, as in the case of mapping out the parameter 
space with MCMC. 



The main reason for using this approximate method 
is speed 



Obviously, the Fisher matrix method has 
its limitations, as for example, was recently discussed in 
(35] . Still we believe that our usage of the Fisher matrix 
approximation - not to obtain the exact bounds but to 
investigate qualitative issues outlined above - is appro- 
priate. 



By expanding the log of the likelihood function to the 
second order in cosmological parameters Si around the 
maximum, one finds the standard expression: 



I X,Y 



dsi 



Coy' 



dC} 



dsi 



(10) 



Here Cov~"^ is the inverse of the covariance matrix (to 
be defined shortly), and the are the power spectra in 
the temperature and polarization channels. We restrict 
ourselves to the temperature, T, E polarized, and cross 
temperature-polarization C, power spectra. They are the 
only ones relevant for scalar perturbations. 



As mentioned before, we will investigate the role of 
other cosmological parameters using the Fisher matrix 



The elements of the covariance matrix [83|, |8J, |85|, 
which give the errors of the corresponding measurements, 
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(13) 
(14) 



Cci{CTi + w^^Bf), (15) 



CcK^^Bf + Wp'S72). (16) 



In these equations, fsky is the sky coverage, Wj,^ and 
specify the detector noise for temperature and polariza- 
tion respectively, and = e'('+^)^f'/(8in2) jg ^j^g bg^jji 
smearing window function. Here 9b is the full-width, half- 
maximum of the beam in radians, wt and wp are the 
quantities characterizing the detector noise level for tem- 
perature and polarization, respectively. For Planck, we 
take Ob = 7', w^^^^ = 40^K', and Wp^^^ = Se^iK'. 

More precisely, the Fisher information matrix Fij is 
defined through derivatives of the likelihood function 
C{x,p) for data x and model parameters p as 



F 



dp^dp^ 



(17) 



The right-hand side of Eq. (fTT)) is averaged over the 
data a;, weighted by the probability C{x,p) of their ob- 
servation in the fiducial model. Given this definition, 
the Cramer-Rao inequality states that the r.m.s. of the 
best-estimator for a parameter p' cannot be less than 

(F~i)", as discussed, e.g., in |85| . 

Of course, to perform the average in general requires 
mapping out the likelihood function over the parameter 
space and we are back to the problem of computing the 
likelihood at O(IO^) points. We do not perform such 
mapping here. Instead, we simply assume the likelihood 
is close to gaussian, compute the matrix in Eq. (fTO|) , and 
use {F~^y^ as an estimate of a lu error on the param- 
eter Si. 

We consider the following set of cosmological 
parameters (w^/cj^, w^, ftdc, N^ff , Ncoupied, n^^as 

As, Yhc,t, w^jf^c, Wdc), replacing Nps with 
N^^^ = Nps + Ncoupied ill our list of parameters. 
The first 10 parameters span the parameter space of our 
MCMC analysis of the last Section, while the last three 
are the new additions. We calculate the derivatives 
in Eq. pO|) by symmetric finite differences about the 
best-fit cosmological parameters from the WMAP year 
three data. We compute the resulting la marginalized 
error on Ncoupied and N'^f^ with all 13 parameters and 



then drop the last three parameters in turn. We also 
consider the effect of fixing the helium abundance Yhc- 

The results are shown in Table IIVI We see that keep- 
ing the curvature and the neutrino mass fixed to zero is 
completely justified. Moreover, the effect of varying the 
dark energy equation of state is also quite small (~ 7% 
for Ncoupied and ~ 9% for N^^^), as is the effect of vary- 
ing Fhc (~ 4% for Ncoupied and ~ 9% for N^H). Thus, 
practically speaking, we are justified in our choice of the 
10 cosmological parameters used in our MCMC calcula- 
tions, as the other three parameters do not change our 
qualitative conclusions about Planck's sensitivity to the 
number of coupled neutrinos. In addition, though the 
Fisher analysis shows that the effects of Irc on the errors 
of Ncoupied and N^^^ are fairly minimal, we include Ifjc 
as a parameter in our MCMC with priors consistent with 
astrophysical observations as described in Sect. IIV A 31 



B. The compensation mechanism and the nature of 
Planck's sensitivity 

We now wish to establish the nature of Planck's sensi- 
tivity to neutrino self-coupling. By this we mean showing 
which effects of self-coupling cannot be entirely compen- 
sated by adjusting the cosmological parameters and how 
the residual differences compare to the corresponding er- 
rors. The latter, as we saw in the previous subsection, 
are set by a combination of the cosmic variance and the 
resolution/sensitivity of the apparatus. Equipped with 
our numerical results, we are now able to explicitly see 
the compensation mechanism in action. 

For this, we turn to Table [1111 The best-fit parameters 
listed in the Table are precisely those which compensate 
the effects of neutrino coupling most efficiently within 
each scenario. Indeed, we recall that while the mock data 
is always generated under the assumptions of three freely 
streaming neutrinos, for each of the scenarios listed in 
the Table, the MCMC code attempts to find the best fit 
possible within the framework of that particular scenario. 
The code varies the values of the cosmological parameters 
until the closest agreement is found. 

The compensation turns out to be very efficient. One 
way to illustrate it is to plot the relative differences of 
the temperature and polarization power spectra in the 
scenarios of coupled vs. freely streaming neutrinos that 
remain after the cosmological parameters are adjusted. 
This is done in Fig. 1111 where we show the situations for 
1,2, and 3 coupled neutrinos. We can see that, unlike the 
differences shown in Fig. [T] which were in tens of percent, 
these residual differences are only at the level of 1-3%. 

We further observe that below the multipole number 
of I ~ 1300 the residuals have no clear structure, while 
above they develop a shape of periodic oscillations, as a 
consequence of the phase shift, Eq. ([7]). 

In terms of the sensitivity of Planck, the physically 
relevant quantities are not the relative differences in the 
C/'s just shown, but the ratios of the differences to the 
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Vary all 
13 params. 


Fix 

uov = Q 


Fix uj^ = 0, 

^curv — 


Fix UOv — 0, flcurv — 0, 
Wdc = -1 


Fix uj^ — 0, ^curv — 0, 
Wdc = -1, Yho = 0.24 


Error on N^oupied 


0.31 


0.31 


0.31 


0.29 


0.28 


Error on N^^^ 


0.38 


0.38 


0.38 


0.35 


0.32 



TABLE IV: The Fisher matrix estimate of the la error on Ncoupied and N^'. The effects of fixing neutrino mass, curvature, 
dark energy equation of state, and the Helium fraction are considered. 
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FIG. 11: The relative residual differences in temperature and 
polarization power spectra for 1, 2, and 3 coupled neutrinos 
(in order of deviation from zero) that remain after the cos- 
mological parameters are adjusted in each of the scenarios (as 
shown in Table llll)l . 



corresponding errors. The errors are given by the ele- 
ments of the covariance matrix, Eqs. (|imi6p . We can 
then consider at each £ the quantities 



iSxj)^ E (ACf)(Cov,)3^V(ACl'), (19) 

X.Y={T,E,C} 

which is nothing but the contribution of a given multi- 
pole £ to the x^- These Figures are very instructive, for 
they show that statistically significant deviations occurs 
for £ < 1300, and have no clear structure. The higher 
multipoles which exhibit clear effects of the phase shift 
in Fig. [TT] turn out to be of relatively minor importance. 

This is further seen when we consider the effects of the 
expected detector resolution and noise on the la errors 
in Ncoupied- In Figs. [131 [TH we show the change in sensi- 
tivity as a function of detector characteristics, obtained 
a Fisher analysis. The rate at which increased resolution 
and lower noise improve the sensitivity to ANcoupied is 
rather modest. Decreasing the noise by an order of mag- 
nitude, for example, improves the sensitivity by less than 
a factor of two. In addition, we consider the sensitivity of 
a hypothetical experiment which measures all multipoles 
with £ < £max within cosmic variance; we see in Fig. [15] 
that sensitivity to multipoles £max S 1000 is sufficient to 
measure ANcoupied = 0.3, and implies that Planck's high 
sensitivity to neutrino free-streaming does not depend on 
the high £ multipoles. It also implies that our qualitative 
predictions for Planck have a high degree of robustness, 
as the sensitivity does not depend strongly on the high £ 
properties of the detector. 

One characteristic of Planck that proves crucial is the 
simultaneous measurement of both temperature and po- 
larization. This makes it possible for Planck to constrain 
Ncoupied down to 0.5 at la level even if Zmax is only 800. 
In comparison, we saw that WMAP3 places virtually no 
constraint on the number of self-coupled neutrinos. The 
difference is the relatively poor quality of the polariza- 
tion measurements at WMAP, as well as the temperature 
measurements for 500 < I < 800, where Planck is cosmic 
variance limited. 



(18) 



XX 



VI. CONSTRAINTS ON NEUTRINO 
INTERACTIONS 



where X runs over T, E, and C, and AC^ are again the 
differences between the power spectra of the reference 
model {Nps — N^-^^ — 3) and the best-fit power spectra 
in a scenario with self-coupled neutrinos. 

The quantities (Sxe)''^ are plotted in Fig. [T^ (the top 
three panels). In the last panel of the figure, we plot the 



Since we have shown that even a single interacting neu- 
trino or a single extra neutrino will be observed or ex- 
cluded to high precision with Planck, we consider the 
implications for models of neutrino-scalar interactions, 
including neutrino dark energy. Before turning to the 
specific models, we first review general constraints on 
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FIG. 12: The quantities (Sxe)^ and (Sxj) defined in Egs.lTSl 
and ll9l They illustrate the ratios of the residual discrepancies 
at each £ to the corresponding errors, as explained in the text. 
The residual discrepancies at ^ < 1300 play a key role. 
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FIG. 13: la sensitivity of Planck (calculated with a Fisher 
matrix) to ANcoupied as a function of the detector noise. A 
noise factor of one corresponds to the sensitivity chosen in the 



main analysis {w. 
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FIG. 14: la sensitivity of Planck (calculated with a Fisher 
matrix) to ANcoupied as a function of the detector resolution. 
The resolution chosen in the main analysis is 7'. 



the couplings g from Majorana neutrino-scalar interac- 
tions which will be possible with Planck data. 

The loss of neutrino free-streaming and production 
of additional relativistic degrees of freedom as a re- 
sult of neutrino-scalar coupling have been considered in 
Ref. [26j . We summarize and reformulate the constraints 
here; a more complete treatment is left to the appendix. 

For a Majorana mass term, the relevant interactions 
^ vv <-> vv, and (if > 2my) vv <-> 0. The 
last interaction is often most constraining on the coupling 
g, and its rate is 



5^ 

16^—"^^ 



(20) 



for < T . This rate can also be described from on 
resonant s-channel scattering; see the appendix for de- 
tails. The m^/T factor accounts for the boost from the 
rest frame of the scalar. This rate increases as the tem- 
perature decreases, and will cause to decay to neutri- 
nos, increasing N^^^ ^ once drops below T, provided 
g> 10-i3(T,,,/m^). 

If, on the other hand, remains relativistic through 
recombination, this same process may tightly couple the 
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FIG. 15: la sensitivity of a hypothetical CMB experiment 
to ANco upied 'is a function of the maximum muhipole, Imax^ 
measured to within cosmic variance. Multipoles with I > Imax 
are discarded. 

neutrino to the scalar, removing neutrino free-streaming. 
The 2—^1 process can only occur in a small region of 
phase space which depends on the angle between the in- 
coming neutrinos, 9 ~ m^/T. In order to isotropize the 
neutrino momentum, therefore, this process must occur 
N = {T/mjiY times. We require then V{vv ^ (f) > 
NH{Trec) I2y|- The equilibration occurs before recom- 
bination ii g ^ 10^^'^ for ~ Tree, and becoming less 
restrictive with dropping by the factor (r^ec/w^)^. 

The exception is if scattering of 0's through a 0^ in- 
teraction may be fast enough to isotropize the scalar mo- 
mentum in between decays and inverse decays. Then 
the constraint is transmuted to T{(f)(j) <-> > T(h'h' 
0) > H{Trec), which implies that if X^T > T^/Mpi, 
where is the coupling of the cj)'^ interaction, and if 
g'^m'^/T > T^/Mp/, then the neutrinos will be strongly 
coupled with the scalars. The equilibration then occurs 
before recombination if 

g > lO-'\Trec/m^) (21) 



[281 1 ; for this process, is an evolving function of red- 

shift as the neutrinos annihilate. We do not consider this 
scenario further in this paper. 

We briefly discuss the constraints from BBN. The pre- 
cise upper bound on the effective number of neutrinos is 
somewhat controversial and dependent on the particular 
set of constraints from data utilized, but is on the order 
of N^ff < 4.5 [l|]. Any additional scalar degrees of free- 
dom brought into equilibrium before BBN will contribute 
SN^ff = 4/7. One scalar alone brought into thermal 
equilibrium before BBN, then, will not upset the current 
bounds. If the theory contains multiple scalars, however, 
there will be constraints already from BBN on the cou- 
pling g from vv ^ (j) and vv <-> 00 coming from demand- 
ing that these processes not recouple until T < Tbbn- 

If these processes are brought into equilibrium after 
the neutrino processes decouple from the heat bath at 
T ~ 1 MeV, they may still affect the value of N^^^f 
measured by the CMB [2^ [23|. For example, it may 
be possible that new particle states {e.g., scalars, sterile 
neutrinos) are populated by recoupling after BBN and 
then annihilate when the temperature drops below their 
masses (but still before the CMB decoupling). In this 
scenario, while the annihilation step occurs at constant 
entropy, the recoupling does not. In the approximation 
when the former occurs at constant energy, the value of 
N^ff measured by the CMB is given by Eq. (8) of [27i] . 
In general, for a given model, accurate analysis of recou- 
pling may be required. 

We summarize the future constraints from Planck 
given no deviation from standard relativistic energy den- 
sity and neutrino free-streaming in Table IVl 



VII. SOME IMPLICATIONS FOR MODELS 
A. Models of Neutrino Dark Energy 



and 



> 10 



-14 



(22) 



2 <-^- 2 interactions offer complementary constraints to 
the decay-inverse decay processes. If both neutrinos and 
scalars are relativistic at decoupling, the rate for lyiy ^ vv 
is given by 



T(vu w) ~ — — r, 



(23) 



and is competitive with vv <-> 00, so that if <? > 10~^ 
neutrinos and scalars will be tightly coupled. If > 
Tree, however, the constraints on g are relaxed by T/rrifj,, 
due to propagator suppression from the scalar. 

We have assumed everywhere that <C Tree- If 
g > 10~^ and > Tree, ^ti^i the neutrinos will anni- 
hilate to the scalars, producing a "neutrinoless" universe 



Models of neutrino dark energy are typically character- 
ized by a coupling between a singlet neutrino and a light 
scalar field of the form given in Eq. [2l An important 
difference between the standard see-saw and the MaVaN 
scenario is that m^v is a continuously evolving parame- 
ter, since the vev of the singlet scalar, (0), varies as the 
universe cools. 

The evolution of (0) results from finite temperature 
effects in the scalar potential; the background neutrino 
energy density acts as a source for the scalar field, and 
may displace (0) from the value dictated by its zero- 
temperature potential, Vb(0): 

1/(0) = Fo(0)+p.(0) 

1 n ,n 77r2 , ml((j))T^ , , 

2 120 24 ' ^ ^ 

where is the light mass eigenstate resulting when the 
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c oust r cLi nt 




uv < — > (j) 


-tree I'f'tp 
Tree ^ f^<p 


g ^ i\j [I ev /rrifp) 
g < 10-" (m^/1 eV)'/2 


R,6iiiovG frGG-strG&ming 
Increase N^-^-^ 


vv ^ uu 


Tree ^ Tn(j> 
Tree '^^(f) 


9 < 10- ' 
9 < 10"' (m^/1 eV) 


Remove free-streaming 
Remove free-streaming 


VU <j)(j) 


< Tree < rUi, 


9 < 


Increase N^-^^ 



TABLE V: Constraints on the effective coupling jr of a Majorana neutrino with a scalar, if the Planck satellite observes no 
extra, non-standard relativistic degrees of freedom or strong neutrino-scalar interactions. 



sterile neutrino is integrated out: 

ml 
= — — 

TUN 

m|_ 



(25) 



This is a self-consistent formalism provided that the ster- 
ile neutrino remains thermally unpopulated (and hence 
effectively integrated out) in the early universe; it was 
shown in [sHi how this can be done with Planck sup- 
pressed operators between the scalar and electrons. In 
addition, in a theory with the full flavor structure of 
three generations, rriy must correspond to the lightest 
mass eigenstate, which is still relativistic in the universe 
today, so that the neutrino dark energy does not clump 
and become unstable [s^l- 

As was shown in Sect. |lTl the mass eigenstate couples 
to the scalar through the Majoran interaction, C = g4>i'i', 
where g ~ \{mD /tumY ■ Minimizing eqn. 1241 



= mo 




(26) 



Hence, we have a temperature dependent coupling g. 

How can MaVaNs be detected with Planck? Planck 
will be sensitive to couplings g, through lyiy <-> vv, of size 
g > 10~^; these couplings may be small enough to lead to 
possible detection or interesting constraints on MaVaNs 
with Planck. We conclude, using mpf from eqn. [21] 



9 ■ 



(27) 



implying that the neutrinos will form a tightly coupled 
fluid down to recombination temperatures, T ~ 0.1 eV, 
unless 



m4. < 10"® eV. 



(28) 



For scalars in these mass ranges, Planck will be able to 
detect MaVaNs. 

What are natural parameters and typical scalar masses 
for the theory? We consider a generic class of MaVaN 
models laid out in [U, [2^ . All energy scales in MaVaN 
models are characterized by the dark energy scale, A ~ 
lO"^'^ eV. All parameters in the model, including m^i 
and UV cutoff of the neutrino sector, are characterized 



by meV scale physics, in contrast to the traditional see- 
saw, which is characterized by a sterile neutrino at the 
GUT scale. Radiative corrections will drive the scalar 
mass to the cutoff of the theory, implying that without 
an unnatural fine-tuning, the scalar mass should not be 
much lighter than the cutoff of the theory at a meV. A 
superpotential which radiatively controls the scalar mass 
is simply the see-saw, eqn. [5J with all the fields promoted 
to superfields: 



W = X(j)NN + yHLN. 



(29) 



Assuming (|A^|) — 0, the superpotential generates no 
tree-level mass term for (j). A mass is generated radia- 
tively for 0, however, through a 1-loop diagram with an 
active and sterile sneutrino in the loop so that 



167r2 



■logAVwD 



(30) 



where A is the UV cutoff of the neutrino sector, typi- 
cally around an eV. We can see that loop corrections will 
generically generate a mass for the scalar, ^ Am^i. 

Besides the loop corrections of Eq. ([50)1 one has to 
keep in mind corrections mediated by gravity. The nat- 
ural size of gravity-mediated corrections to the scalar 
mass from supersymmetry broken at the TeV scale is 
^ F/Mpi ~ lO"'^ eV, or some five orders of magnitude 
above the bound in Eq. 

A and mD are themselves constrained by the require- 
ment that MaVaNs produce dark energy. In particular, it 
was shown [2^ that a mass term appears for |7V1 through 
radiative corrections of size — A^m|,, as it does for \cf>\, but 
with the opposite sign. This implies a false minimum at 
(|A^|) = 0, with a true minimum at (|iV|) tod/A ; 
the energy difference between the true minimum of the 
potential and (|A^|) = is the dark energy density, 



PA 



^i,/A^ 



10-1" eV^ 



(31) 



and, using Xmo from naturalness, we find that 

the only underlying free parameter is the coupling A: 
^ XniD A'^/^IQ-^'^ eV. In terms of the Lagrangian 
parameter A, then, Planck will be able to constrain the 
coupling to be 



A < 10 



'5.5 



(32) 



Since A is a parameter typically envisioned to be not 
much smaller than 1, the reach of Planck is considerable. 



21 



Of course, these constraints on Lagrangian parameters 
depend on naturalness requirements; the direct (more 
model independent) constraint is expressed in terms of 
the scalar mass in eqn. 

The same parameter A enters into matter effects in 
neutrino oscillation experiments, studied in, e.g., [ssLlsol. 
S im, nil. There it was shown that the shift of the 
neutrino mass due to earth size matter effects is 

(33) 

where oat is a Planck suppressed coupling between nu- 
cleons and 4>, and p^v ~ 3g/cm'^ is the earth energy 
density. Neutrino oscillation experiments in earth typ- 
ically cannot reach sensitivities to mass scales smaller 
than 10~^ — 10~^ eV. From this, we can can see that 
Planck and neutrino oscillation experiments are compli- 
mentary probes of neutrino dark energy. Neutrino os- 
cillation experiments are efficient at detecting effects in 
the low 7710 region of parameter space; using the natu- 
ralness requirement, ~ Amu, and the dark energy 
requirement, mjy/X^ ^ p\, 

— (^)fe)ra"' 

Planck, on the other hand, can probe the large region 
according to eqn. [28l 

B. Late time neutrino masses 

In models of late time neutrino masses [1^ [13] , (p car- 
ries a C/(l) charge, which has a low scale / of symmetry 
breaking. Since the mass eigenstate in the operator of 
the form gcfivv only acquires a mass once the U{1) is bro- 
ken (exactly as in the Higgs mechanism), the neutrinos 
only gain their masses at late time. For a low-energy 
Majorana mass term, we expect g ttIu/ f , as shown in 
Sect, ini We conclude then that Planck will be able to 
probe scales of neutrino mass generation: 

^1 MeV < 2m^ , . 

^ - 1 TeV > 2m^ ' ^ ' 

assuming '-^0.1 eV. Hence Planck may probe TeV 
scale neutrino mass generation. 



The current data disfavors 3 non-free-streaming neu- 
trinos at the 3.5cr level, however, this exclusion is not 
without some controversy. The exclusion comes from the 
SDSS Main data set (and also from the Lyman-a for- 
est data), which favors five free-streaming neutrinos. We 
have seen that the LRG SDSS dataset instead favors close 
to three free-streaming neutrinos, in agreement with the 
standard model expectations, but has little sensitivity to 
coupled neutrinos. 

We have shown that Planck will be able to resolve the 
present controversy. Planck alone will be capable of ex- 
cluding a single interacting neutrino at the 4.2cr level. 
Thus we expect sensitivity to models of neutrino mass 
and neutrino dark energy which is currently unavailable. 
Models of Mass- Varying Neutrinos, for example, in many 
regions of parameter space, predict only a single strongly 
coupled neutrino. Likewise, models of neutrino mass may 
have a fiavor structure such that only one or two of the 
neutrinos is strongly enough coupled to the scalar to re- 
move neutrino free-streaming at the CMB epoch. For the 
MaVaNs models, the reach of Planck will be considerable, 
probing the scalar masses down to 10~^ eV, thus 

covering the parameter values favored by naturalness. 

In the case of extra thermalized neutrinos or scalars, 
the current constraints depend on which data sets are 
chosen in the analysis, and in any case the errors en- 
compass multiple neutrinos; as a result, the current data 
are not conclusive on the total radiation density at the 
CMB epoch. By contrast, Planck alone will determine 
the radiation content to AA'^j^ — -o 3- 

The nature of the neutrino mass generation mechanism 
and beyond-the-standard-model neutrino interactions is 
currently unknown. The coupling between a neutrino 
and a light scalar in a Majoron interaction term of the 
form g4)i>v can be as small a& g ^ 10"^'^ (depending on 
the mass of the scalar) while still being large enough to 
remove neutrino free-streaming and make an imprint on 
the CMB. The constraints (and potential for future de- 
tection) considered in this paper should be a guide to 
model builders when considering potential signals from 
future CMB experiments. It has been challenging to test 
neutrino properties in terrestrial experiments, and the 
nature of the mechanism which generates neutrino mass 
remains unknown; future CMB experiments may have 
the tools to unmask its origin. 



VIII. CONCLUSIONS 

We have considered in this paper the effect of neu- 
trino interactions on the CMB, both through removal of 
neutrino free-streaming and the thermalization of extra 
relativistic degrees of freedom. The constraints from the 
current data are shown in Figs. SI [S] and [51 those expected 
to be obtained with Planck are depicted in Figs. [5] and 
ini These constraints are summarized in Fig. 1161 
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FIG. 16: Summary of the bounds on the numbers of free-streaming and coupled neutrinos from the current data (left) and 
the simulated Planck data (right). For the case of the current data, the Large Red Galaxy (LRG) dataset from SDSS yields 
a region that differs significantly from what is obtained with the Main SDSS dataset. The combined fit is also shown. The 
sensitivity of Planck to the numbers of both freely streaming and coupled neutrinos will be a dramatic improvement over that 
of all present-day experiments combined (compare "Planck only" and "All current" on the right). If Planck's data is combined 
with today's data, further improvement on the constraints are expected. 



APPENDIX: ON RESONANT PRODUCTION 

In this section we perform a more rigorous derivation 
of the results of Sect. VI. 

Consider a "test" neutrino traveling through a neu- 
trino gas. We need to know how far it travels before ap- 
preciably changing the direction of its momentum, i.e., 
given the initial momentum of order T, when the parti- 
cle's momentum in the transverse direction becomes or- 
der T. The rate for this process needs to be compared to 
the expansion rate of the universe, given by ~ T'^ /Mpi in 
the era of radiation domination. 

Let us consider the scattering of two neutrinos, vv — > 
viy, and start with the t-channel exchange of (f>. In this 
case, for light (j) (w^ ^ T) the cross-section is strongly 
forward peaked (the well-known property of Ruther- 
ford scattering). The situation is well-known in plasma 
physics. The build-up of transverse momentum happens 
mostly as a result of many small-angle scattering events, 
rather than a single large-angle event. The multiple scat- 
tering events result in a random walk process for the 
transverse momentum px, so that the rate of change of 
is given by 



dpj 
dt 



d{cos 0) 



da 

dcosi 



-nVrelPe- 



(A.l) 



Taking Vrei ~ 1, T'^O'^, n and da /d cos 6 

<7V(167r)r2(T26'2 



mi) ^, we get 



(A.2) 



The small-angle behavior of the cross section is regulated 
by the mass of the scalar field. In plasma physics, the 



logarithm similar to that in the above equation ("the 
Coulomb logarithm") is regulated by plasma screening 
effects. 

The rate for the process needs to be compared to x 
T'^/Mpi ("momentum exchange of order T by the time 
when the temperature equals T"). We get that neutrinos 
are streaming freely at recombination when 

5 < (Trcc/Mpi)'/^ X (167r/ ln(r,ec/m^))'/^ (A.3) 

The 1 /4 power makes the coefficient (the second term) of 
order one for a wide range of m^, hence for the purpose of 
this estimate we can drop it. The order of magnitude is 
set by the ratio (Trec/Mpi). For the temperature 0{eV) 
one thus finds the neutrinos are not coupled by the t- 
channel scalar exchange if 



(A.4) 



The 1/4 power also means that when the t channel pro- 
cess is dominant it is reasonable to treat neutrinos as 
tightly coupled on all scales, not just the sound horizon 
at the CMB decoupling. Indeed, even considering scales 
100 times smaller than the sound horizon only changes 
the bound on g by about a factor of 3. 

If TO0 > T, the situation is different. The momentum 
exchange occurs via large-angle scattering and the cross 
section is set by T^/ml, so that 



5> (Trcc/Mpi)'^^("^0/?lcc) 



(A.5) 



The interaction decouples as the temperature drops, a 
characteristic feature of nonrenormalizable interactions. 
Indeed, at temperature T <C the (p field can be inte- 
grated out, resulting in an effective 4-fermion vertex. 

We now consider vv — > vv scattering via the s-channel 
process. In this channel, an important possibility is the 
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scattering on-resonance. When it is open, it can increase 
the sensitivity to g by orders of magnitude, as we will see 
next. 

For 7710 < r, the resonance condition is satisfied when 
the relative angle 9i2 between the two incoming neutri- 
nos, 1 and 2, is small, so that s = i?i£^2(l— cos6'i2) — m^, 
or 012 ~ m^/T. The relative velocity v^ei is given by 9i2. 
The cross section in the center-of-mass frame is given by 
the relativistic Breit-Wigner form. 



(A.6) 



where T is the width of the resonance, given by the rate 
of decay of (j> in its rest frame, — g^m0/(167r). The 
cross section has a peak value ^ Tr/m'^. In the frame of 
the universe, the cross section is reduced by a factor of 
m^/T = 9i2 (since it is boosted in the direction perpen- 
dicular to the direction of the colliding momenta in the 
center-of-mass frame). 

Collecting all the factors, we obtain 



712"-resC'12 



T^9 



12i 



(A.7) 



where rires is the number density of neutrinos on which 
the test neutrino scatters resonantly. It can be estimated 
from the phase space as ~ Tm^T^ = g^Tm?^/ {IQn). 
Hence, 



dt 



167r T 



(A., 



The result in Eq. (jA.SP has a simple physical interpre- 
tation [2^, as a resonant production of the (f> particles, 
— s- (f), completely analogous to the Z-boson produc- 
tion at LEP. Indeed, the rate can be most easily com- 
puted in the opposite direction: the decay of (f) boosted 
into the frame of the universe is g^rn^/T/ {16tt). To get 
the momentum exchange rate, one recalls that for each 
scattering event, one has pt w^. 

The crucial feature of the resonant process is the de- 
pendence on g^, rather than g*^ since as we saw the 



reference values of g are in the 10 ^ range. Instead of 
Eq. (jA.3|) one has 



9 < (4^^/')(r,ec/MpO'/' X (T,ec/m0)2 



(A.9) 



For TO0 > 10^"^ eV this constraint is stronger than that in 
Eq. (|A.3p . Otherwise, the integral over the nonresonant 
part of the phase space dominates, giving a bound similar 
to Eq. (|X3|l . 

So far, our discussion has been in complete agreement 
with Ref. [2^. We now point out some differences and 
make additional observations. 

As noted in [2^, in the models where the light scalar 
(j) is the Goldstone boson there is a possibiHty to produce 
(j) via the s-channel diagram with the heavy scalar p (the 
radial mode) in the intermediate state, p —> cfxj). 

Just like the process vv ^ <f> ^ vu considered above, this 
process is most efficient at producing the phi particles on 
resonance, i.e., at the point in the history of the universe 
when the temperature equals ~ rup. Unlike the result 
given in Eqs. (28) and (29) of [1^ we estimate the rate 
of this process treating it as I'z^ ^ p 



9 

327r' 



(A.IO) 



Comparing it with T"^ /Mpi — m? JMpi we find that unless 



g < 4(2^)1/2 (m,/AfpO'/' 

< io-"(mp/r,,e)^/' 



(A.11) 



the Goldstone degree of freedom will be populated. If 
this occurs before the neutrinos decouple (i.e., if rUp ^ 1 
MeV), the number of effective neutrino degrees of free- 
dom at BBN will be 3 -I- A/7Ng where Nq is the number 
of Goldstone fields. For a single Goldstone this may be 
still acceptable, but more than one would cause tension 
with the BBN bounds on N^, as well as with the projected 
bounds of Planck, as we have shown here. 
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